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Abstract 

In 1998, Khodkar showed that the minimal critical set in the Latin 
square corresponding to the elementary abelian 2-group of order 16 is 
of size at most 124. Since the paper was published, improved methods 
for solving integer programming problems have been developed. Here we 
give an example of a critical set of size 121 in this Latin square, found 
through such methods. We also give a new upper bound on the size of 
critical sets of minimal size for the elementary abelian 2-group of order 
2": 4" - 3" 4 - 2" - 2""^ We speculate about possible lower bounds for 
this value, given some other results for the elementary abelian 2-groups of 
orders 32 and 64. An example of a critical set of size 29 in the Latin square 
corresponding to the elementary abelian 3-group of order 9 is given, and 
it is shown that any such critical set must be of size at least 24, improving 
the bound of 21 given by Donovan, Cooper, Nott and Seberry. 

1 Definitions 

A Latin square L of order n is an n x n array of entries such that 

each row and column of L contains each of n possible elements exactly once. In 
this paper these n possible elements will be the set {0, . . . , n — 1}, as we will be 
using bit-wise logical operators later. A subsquare S of order m from L is an 
array of m x m entries from L, not necessarily contiguous, such that each row 
and column of S contains each of the m possible elements exactly once. 

A critical set C of L is a subset of L such that L is the only superset of C 
which is a Latin square (is uniquely completable) and no subset of C has this 
property. 

We now review the definition of a Latin trade (see for example jSj). Let P 
and P' be two partial Latin squares of the same order, with the same size and 
shape. Then P are P' are said to be mutually balanced if the entries in each row 
(and column) of P are the same as those in the corresponding row (and column) 
of P'. They are said to be disjoint if no cell in P' contains the same entry as the 
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corresponding cell in P. A Latin trade (also known as a Latin interchange) I is 
a partial Latin square for which there exists another partial Latin square /', of 
the same order, size and shape with the property that I and I' are disjoint and 
mutually balanced. The partial Latin square /' is said to be a disjoint mate of 
I. A Latin trade of size 4 is known as an intercalate. An example of a Latin 
trade is given below. 
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The following lemma clarifies the connection between critical sets and Latin 
trades. 

Lemma 1 A partial Latin square C C L, of size s and order n, is a critical set 
for a Latin square L if and only if the following hold: 

1. C contains an element of every Latin trade that occurs in L; 

2. for each {i, j; k) G C, there exists a Latin trade Ir in L such that IrHC = 

i; k)}. 

Proof. 

1. If C docs not contain an clement from some Latin trade /. where / has 
disjoint mate /', then C is also a partial Latin square of L' = (L \ /) U /'. 
Hence C is not uniquely completable. 

2. If no such Latin trade Ir can be found, then the position {i, j; k) may be 
deleted from C and C \ {{i, j; k)} will still be uniquely completable and 
thus a critical set for L. 



The elementary abelian m-group of order m" is the group multiplied 
with itself n times, denoted in this paper by L{m"'). Elsewhere, it has also been 
denoted by Z^. 



2 Integer Programming 

To find a critical set C in the elementary abelian 2- group of order 2", it thus 
suffices to solve an integer programming problem. If T is the complete set of 
Latin trades in a Latin square L, then the optimal solution to the following 
integer program is the size of the minimal critical set in L. 
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Minimize : Y^xgl^^ 
Subject to : 

For each T G T, ^ > 1 

Throughout this paper, C^; is 1 if x G C and otherwise. 

In Khodkar uses this method to prove that the minimal critical set for 
L(8) is 25 and finds a critical set of size 124 in L{16). Even for L{8), it is 
a difficult task to find all the Latin trades in this square, and so a subset of 
the Latin trades must be used. For L{8) Khodkar chose all the subsquares 
isomorphic to L{2) and L(A) because a linear program similar to the above can 
be written as follows. Ii denotes all the subsquares isomorphic to L{2) in L{8) 
and I2 denotes all the subsquares isomorphic to L{4) in L{8). 

Minimize : E^eL(8) 
Subject to : 

For each Teh, '^C^>1 
xeT 

For each Teh, '^C^>5 

xeT 

As explained by Khodkar, although this does not cover all the Latin trades in 
L(8), this works because the inequalities above make each subsquare isomorphic 
to L{4) in L{8) uniquely completable. Since a critical set of size 25 can be found 
in L{8) and the optimal solution to the above integer programming problem is 
also of size 25, we know that the minimal critical set size in L{8) is also 25. 

Khodkar used the CPLEX solver W, version 6.5, for solving these integer pro- 
gramming problems. For solving integer programs, CPLEX uses a branch- and- 
bound or branch-and-cut approach. Although CPLEX can find provably opti- 
mal solutions, for finding "good" solutions quickly a non-deterministic search 
is much better. We use this approach to find a critical set in L{16) and much 
better solutions to similar integer programs for L(32) and L(64). 

Walser extended the WalkSAT algorithm, a satisfiability (SAT) solver using 
local search |12| . from solving SAT problems to solving integer linear program- 
ming problems. The resulting program is called WSAT(OIP) [2- We consider 
the following problem, also from the Khodkar paper, where is the set of all 
subsquares in L{16) isomorphic to i(2"). 



Minimize : ExeL(i6) 
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Subject to : 

For each T G h, ^ C^; > 1 

For each Teh, ^ C^, > 5 

For each Teh, ^C^>'25 

xeT 

Khodkar found a sohition of size 124 to the above problem with CPLEX 
which was a critical set for L{16). He noted ^ that CPLEX did find smaller 
solutions which were not critical sets. This is the case because the inequalities 
above do not necessarily force the subsquares isomorphic to L{8) to be uniquely 
completable, though they force subsquares isomorphic to £(4) and L{2) to be 
uniquely completable. Also, other Latin trades not specified in the program are 
missing from L(16). 

Using the CPLEX 8.0 solver on this problem, it took several hours on an 
Pentiuni-4 system at 1.5 Ghz, with a variety of options, to find a solution of size 
118. In contrast WSAT 1.105 can find solutions of size 112 to the above problem 
in a few seconds. The disadvantage is that it cannot prove such solutions are 
optimal. We note the right hand side (RHS) of the inequalities for all such 
solutions of size 112 found are as follows: 512 with RHS 1, 208 with RHS 2 or 
3, 32 with RHS 4; 72 with RHS 5 or 7, 40 with RHS 9, 8 with RHS 10, 11 or 
12; 28 with RHS 28; and 16 with RHS 26 or 30. Such symmetry and universal 
divisibility by powers of 2 leads the author to conjecture that 112 is in fact the 
optimal solution to the above problem. Also, in all such solutions of size 112, 
each of the subsquares isomorphic to L{8) is uniquely completable. 

It is simple to find critical sets of size 121 using this method. Here is an 
example: 
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Methods to try to improve this bound such as modifying the RHS bounds 
of 5 and 25 to various other values ranging from 6 or 7 for the I2 trades and 
26 to 30 for the trades have been unsuccessfuL Similarly finding the missing 
trades in solutions of size 25 or more to the third set of inequalities has not 
aided much. 

If constraints are added to the program to ensure that each row and column 
contains either 7 or 8 entries, and each element occurs either 7 or 8 times, a 
critical set of size 121 can still be found. However this is still the best known 
solution. To improve this result, the author suggests that a number of methods 
be tried. 

• The structure of the problem best suits a 64-bit CPU - any subset of 
L(16) can be represented as 4 64-bit words. To reduce the size of a known 
critical set, perhaps the same method as jS] could be used - for increasing 
n, attempt to remove n + 1 entries while adding n entries, and check for 
unique completion afterwards. 

• An approach similar to the paper ^ which concerned Steiner triple sys- 
tems can be used. That is, if we find that our integer program solver 
outputs a solution which is not uniquely completable, we determine which 
trades in -/j(16) are "missing" from the solution, use the automorphism 
group of L{1Q) to generate copies of all such trades, "minimize" this list, 
add them to the integer program, and continue. 

• Recently, Margot |7| wrote about pruning 0-1 integer programs which 
have a high degree of symmetry, and used this method to prove that no 
4-(10,5,l) covering design with less than 51 sets exists. Although the 
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problem for L{16) may have too many variables for an optimal solution 
to be found, there is hope that this method may find a better solution. 
Similarly, the integer programming methods of Applegate, Raines, and 
Sloane JHl may be applicable here. 

Applying the above integer program to L{32) the best solution found so 
far has been of size 546. Using just intercalates, a solution of size 514 has 
been found. It has not been computationally feasible to determine whether or 
not either solution has unique completion. The best bound currently known 
is obtained by a method similar to that of Donovan, Fu and Khodkar [S]; we 
take the critical set of size 4^ — 3^ = 781 as constructed by Stinson and van 
Rees P], add in the main backwards diagonal and the last row and column, and 
recursively remove entries beginning in the top left-hand corner. This results in 
a critical set of size 658. 

Given the evidence above, L(32) seems not to have a critical set of size less 
than or equal to Ghandehari, Hatami, and Mahmoodian ^31 found that 
there exists a Latin square L such that the minimal critical set in L has size 

— (e + o(l))n^/'^. The author conjectures that the smallest order n for which 
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there exists a Latin square L of order n with no critical set of size less than ^ 
is n = 32, in the square L(32). 

Similarly, the best solution to the above integer program for L(64) has been 
of size 2470. 

With these results in mind, we conjecture the following bound. 

Conjecture 1 If C is a critical set of minimal size size in _L(2"), n > A, 
121.4"--* < |5| < 5"-i. 

3 A new bound for the size of a minimal critical 

set in L(2") 

We can find a new general bound on the size of critical sets in i(2") by taking 
the same Stinson and van Rees construction for a critical set in i(2"), adding 
3 entries and deleting 2" + 2"~^ — 1 entries. This results in another critical set. 

Theorem 1 If C is a critical set of minimal size in i(2"), |C| < 4" — 3" + 4 — 
2" - 2"-^ 

Proof. We describe a construction isomorphic to the one given by Stinson and 
van Rees for a critical set C in L{2"). Formally, if we consider the critical set 
C in L(2"), the entry at {i,j) G C if and only if {i & j) ^ 0, where & is the 
bitwise logical AND operator. Since we are dealing with i(2"), the element at 
(i, j) is i A j where A is the logical XOR operator. 

Thus C = {ii,j]iA j)\{ikj) ^ 0}. Consider D = {C U {(0,0;0), (0,2" - 
1; 2"-l), (2"- 1,0; 2"-l)})\({(2"-l, x; 2"-l-a;), {x, 2"^1; 2"-l-a;)|2"-i < 
a; < 2" - 1} U {{x,x; 0)|2"-2 <x< 2"-* - 1}). 
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D has unique completion as we first add back in the entries {{x, x; 0)|2" ^ < 
X < 2"-i - 1}, followed by the entry (2" - 1, 2" - 1; 0), then the entries {(2" - 
l,x;2" -1-x), (a;,2"-l;2"-l-x)|2"-i < a; < 2" -1} in any order, resulting 
in a superset of C. 

To prove that every entry in C is necessary, we say that an intercalate 
/ C L{2") proves the necessity of an entry x & C ii I CiC = {x}. If we can find 
such trades for every entry in C, then C is a critical set, by Lemma 1. 

We prove the necessity of the two additional entries (0, 2" — 1; 2" — 1), (2" — 
1, 0; 2" - 1) by considering the two intercalates: {(0, 2" - 1; 2" - 1), (0, 2"-^ - 
l;2"-i-l),(2"-i,2"-l,2"-i-l),(2"-i,2"-i-l;2"-l)} and {(2" - 1, 0; 2"- 
1),(2" - l,2"-i,2"-i - l),(2"-i - l,0;2"-i - l),(2"-i - 1, 2"-i, 2"-^, 2" - 
1)}. For the entries {{x,x;0)\0 < x < 2"^^ — 1}, the following intercalate 
proves the necessity of each entry: {(x, x; 0), (x, 2"^^; x + 2""^), (2"~^, x; x + 
2"~^), (2"~^, 2"~^; 0)}. For any other entry x £ D, the intercalate which suffices 
to prove the necessity of x in C also suffices for D. Therefore D is a critical set 
of the required size. ■ 

4 The elementary abelian 3-group of order 9 

We create an integer program IPi similar to the above by adding constraints 
based on all the Latin trades on 6 or less rows, columns, or elements which have 
size less than or equal to 20. We found, using WSAT, the following critical set 
C29 of size 29 in the elementary abelian 3-group of order 9, L{9). Donovan, 
Cooper, Nott and Seberry JI] had found that any critical set L{9) contained 
at least 21 entries. 
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We create a new integer program IP2 with all 324 Latin trades of size 6 in 
the Latin square, and add symmetry considerations as follows. We call rows 
to 2, 3 to 5, and 6 to 8 of L(9) R2 and R3 respectively, and similarly call 
columns to 2, 3 to 5, and 6 to 8 Ci, C2 and C3. Ei, E2, and E3 represent 
elements to 2, 3 to 5, and 6 to 8. There are nine subsquares defined by the 
row and column sets. If C is a critical set in L(9) we refer to these subsquares 
as C{Ri,Cj). 
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Theorem 2 If C is a critical set in L{9) then there exists another critical set 
D mL(9) such that \D{Ri,C\)\ > \D{R,,Cj)\ for all (ij) e {1, 2, 3} x {1, 2, 3}. 

Proof. In C, wc can permute Ri with R2 or i?^, and Ci with C2 or C3 until 
the square defined by i?i and Ci contains at least as many entries as any of the 
other eight defined subsquares. Call this property RC. 

Each of these permutations is based on the fact that for any critical set 
C in L(9) there exists an isomorphic critical set in L{9) after any of these 
permutations are performed. This is true because after the permutations above, 
we can permute R2 with i?3, C2 with C3, or Ei, E2 and £"3 to produce a critical 
set isomorphic to C in i(9). For example, for the critical set C29 above, we can 
swap Rx with i?3, then we can swap Ei with E^ and C2 with C3 to produce a 
critical set in L(9) with the property RC. If for another set, we swapped C\ 
with C3, we could then swap R2 with R^ and Ex with to obtain a critical 
set in i(9). Similar arguments apply for other permutations of rows and/or 
columns. ■ 

Thus IP2 has another eight constraints for the property RC , as in Theorem 
2. Solving IP2 with CPLEX, we find that any critical set in this square is of 
size at least 24, improving the bound of Donovan et al. 

Theorem 3 If C is a critical set of minimal size in i(3"), n > 2, \C\ > 
24.9"-2. 

The inability of WSAT to find a solution of size less than 26 strongly suggests 
that 26 is the size of the smallest solution to IP2 ■ 

Solutions of size 28 to the integer program IPi above have been found; 
some such solutions have trades with more than 6 rows, columns, and elements 
missing. Thus a better approach, perhaps based on one of the three suggestions 
above, is needed to prove that 29 is the optimal solution, or that a better 
solution can be found. It would be symmetrically aesthetic if a critical set of 
size 27 could be found with exactly 3 entries from each of the 36 3 x 3 subsquares 
here. Over 540,000 solutions can be found for an integer program with precisely 
these constraints (constraints based on the 324 trades of size 6, plus constraints 
ensuring there are exactly three entries in each subsquare). 

Conjecture 2 The critical set of minimal size in L{9) is of size 29. 
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